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ERRATA IN WIS-TCI-469 


Location 
line 3 


Bottom 


line 5 
Bottom line 
line 2 
Ref. 19(a) 
Ref. 36 


Change to read 
o 

After 9.4 A), add sentence: 

20b 

(However, there is a stated uncertainty of 

o 

some 25% in the Buck-Pauly potential at 9 A.) 

Add paragraph: 

Finally, the precision attainable in the 

2 

determined parameters (such as e r m and £ ) 
is limited by the degree of precision of A^, ..., 
A^, from Table II. For example, using Model I 
and Eq. (17b), the statistical uncertainty in 
the determined value of £ will be at least + 8% 
under the best of conditions. 

Degree Fitted^ 

°2 

case E , 8 = 1 A 

what similar in shape, ... 

J. Chem. Phys. J>6 (1972), in press. 

(a) U. Buck (etc.) (b) U. Buck, private communi- 
cation (1972); to be published. 
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ABSTRACT 

With the advent of improved experimental data on absolute total cross 

sections Q(v) for atom-atom collisions and their velocity dependence, on 

the glory undulations and the transition to high-velocity behavior, it 

is timely to reconsider the problem of inversion of such data to yield 

information on the interatomic . potential. In the absence of additional 

data in the form of differential cross sections there is a limit to the 

amount of information available from Q(v) even when observations of good 

accuracy ( e.g . , ± 0.257.) are in hand over an extended energy range ( e . g . , 

3 

from "thermal" energies upward by a factor of _> 10 in relative kinetic 

energy). A number of commonly used procedures for data inversion are no 

longer adequate to deal optimally with the high quality experimental 

results now becoming available. The present paper attempts to develop 

improved methods for data utilization/ which take full advantage of the 
accuracy of the experimental Q(v) measurements. 
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I . INTRODUCTION 

It was recognized in 1934 by Massey and Mohr^ that total cross sections 
for the elastic scattering of atoms by atoms could yield information on 
the interatomic forces, particularly the long-range part of' the inter- 
action potential. At about the same time the experimental (beam) techniques 

for making such measurements (as originated by Born,^ a Bielz^ and Knauer 
2 c 

and Stern in Germany during the previous decade) were brought to fruition 

3 

quantitatively by Mais, Rosin and Rabi in the U.S.A., by Sasaki, Nishibori, 

4 5 

Kodera and co-workers in Japan and by Fraser and Broadway in England. 

In 1935 the status of the subject as a "new kinetic theory of gases" was; 

established, in a summary by Rabi, which complemented an overall review 

6b 

on molecular beam scattering by Guillemin. In 1936 Massey and Bucking- 
ham^ showed how existing total cross section data for alkali-rare gas 
systems could be utilized to ascertain the magnitude of the long-range, 
London dispersion "C^ constant," i.e . , the coefficient of the asymptotic 

g 

attractive potential V(r.) ~ -C g /r , where for S-state atoms s = 6. These 
fexperimentally-derived coefficients compared well with theoretical estimates 
based on the S later-Kirkwood-He llmann (SKH) approximations. 

A rather small number of experimental studies on total cross sections 
ftere carried out over the next two decades, however. One of the more 
extensive of these was a series of measurements (and theoretical cor- 
relations) of total cross sections for the scattering of numerous atom- 

g 

atom and atom-molecule systems, reported in 1959. All of the experiments 
lip to this time had involved the scattering of thermal (Maxwellian) beams 
by thermal "target" gases, so that the resulting cross sections were all 
heavily velocity-averaged. However, the interesting possibility of 
verification of the asymptotic interatomic force law, i.e . , determination 
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of the inverse power s in the long-range potential, from the velocity- 
dependence of the cross section had been known since the appearance of 
the original Massey-Mohr (MM) equation,^ - which can be written in the form 

x 

= fVi n ^ (^sM\r) S ' (i) 

Here the constant of proportionality -|b^(s) is a known, slowly-varying 
function of s , and v is the relative velocity. 

In 1960 the first experimental study of the velocity-dependence of 

9 

the total cross section was reported, by Pauly, for the K-^ system. 

Here s was found to be 6 (within an uncertainty of about 10%) , con- 
firming the theoretical expectation, thus indirectly establishing con- 
fidence in the procedure for determining values from "absolute" values 
of thermally-averaged total elastic scattering cross sections. 

Three sources of doubt remained, however, in the resulting potential 
constants C^. One involved the possibility of systematic experimental 
errors in the cross sections. This suggested itself from the fact that 
many of the experimentally-derived constants were significantly larger 
than theoretically estimated (SKH) values, although ratios of values 
were in excellent agreement with theory. ^ The source of this error^ 

1 la 

was determined and previous experimental results corrected as required. 

The revised constants accorded well with SKH-approximated values, and 
even better with the more rigorous perturbation- theory results of Dalgarno 
and co-workers.^ 

Another source of doubt was the validity of the MM approximation in 
general, and in particular the magnitude of the constant -j?^( s) in Eq.(l). 
Even on the assumption of a pure r ^ potential, several different approxi- 



mation formulas had been derived in the literature, each having employed 

somewhat different assumptions. All of these yielded the same functional 

form of the velocity dependence as the MM equation, but slightly different 

coefficients p(s). Table 1 lists the values of p(6) thus obtained. Not 

14 

included are results of three other approximate treatments, based on 
the application of the uncertainty principle, which lead to the same MM 
form but whose p values are inherently less accurate. 


Table 1. Values of p(6) according to different approximations 

, 1 13a _ , -irj .13b , 13c 

Massey-Mohr Firsov Schif f Landau -Lrf shitz 

p(6): 7.547 8.037 8.083 8.083 


On the basis of a comparison with sample exact calculations (for a pure 

_ 12 3_3d 

repulsive r -power potential) it was concluded that the Schiff- 1 

Landau-Lifshitz (SLL) coefficient p„ TT was to be preferred over p,„,. 

SLL MM 

However, the entire range of p(6) values of Table 1 is only 7.1 °L, so that 
in any case the uncertainty in the derived Cg from a given Q is confined 
to a span of ca. 18%. i 

A third source of doubt was the possible influence of the next higher- 

order terms in the long-range potential, e . g . , the r ^ retardation term 

“8 

and/or the r dipole-quadrupole contribution. A perturbation-like treat- 

15 

ment of these effects on the total cross section has been carried out, ; 

which indicated that the more important effect (at "thermal" conditions) 

-8 

was due the r term, the incremental cross section being approximately 

where p = Cg/Cg . Recent theoretical calculations of P by 
16 

Davison, though significantly larger than values from the older litera- 

16 —8 

ture, have led to the conclusion that the r 


term makes only a small 
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contribution to the cross section at thermal energies for the typical 

systems investigated. Neverthe les s , the presence of this extra contribution 

to the cross section, which is to a first approximation velocity-independent, 153 

-2/5 

could make itself known via a deviation from the v form of Q(v) 

1 7 

(cf. Eq.(l)) This type of deviation was observed by Beck and Loesch 
for the K-Kr and K-Xe systems. The direction and magnitude of the effect 
was in accord with that expected on the basis of Refs. 15a and 16. How- 
ever the results cannot be considered to be definitive evidence for the 
C g contribution, since any additive attractive contribution to the potential 

over and above the asymptotic form V ~ -Cg/r 6 would introduce a quali- 

-2/5 

tatively similar deviation from the v velocity dependence. It is 

therefore difficult to extract from Q(v) data definitive information on 

- fi 

such "correction terms" to the r potential tail. 

18 

Schlier and co-workers have shown by means of calculations based 

on a flexible model potential (termed "realistic") that the Q , averaged 

over the glory undulations, can differ significantly from the Q based 

upon the Cg assumed for the model potential. This is presumably due to 

the fact that the range of r probed by the Q measurements (over the 

velocity range considered) extended inward to smaller separations than 

those for which the potential could be well approximated by its asymptotic 

form, i.e . , in the region of IT near (Q/ZTT) 1 ^ 2 , v(r) ^ -Cg/r 6 

A clear experimental indication of this difficulty is seen from the 

19 

recent work of Pauly and associates. For the systems Na-, K~, and 

2 0 

Cs-Hg for which a full "Buck-inversion" of all scattering data has 

yielded the "true" potentials, 19a > 20b Buck et al . 19a showed that even at 

their largest reduced separation z = r/r = 2 , where Ivl /c ^ 0.015 

m 7 r 1 f > 

“8 

the Cg r and higher terms constitute a significant contribution to the 
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potential. Specifically, for Na-Hg, using the best theoretical estimate 
2 1 l *"6 

for C r one calculates that the C, r term accounts for only slightly 
6 b 

than half the total potential^^ at this separation (of 9.4^). 


more 


22 


This result, if found to be representative, would call into 

question the practical usefulness of the long-range perturbation expansion 

23 

in reciprocal powers of r . Certain and Bruch, reviewing the question 
of the validity of the long-range expansion, indicate that the higher 
multipole terms ( i.e . , Cg, C^q etc.) become important at about the same 
r as the exponential overlap terms. They state that "it rarely makes 
sense to include higher order multipole terms while neglecting overlap 
terms'.'. One might say, that when the leading term does not suffice, the use- 
fulness ’ of the series expansion iis doubtful. Nevertheless one must be careful not to 
"throw out the baby with the bath water". At present, total cross sections 
constitute the only "direct", source of information about the C constants 

and have, on the whole, yielded values in fair accord with theory. 

24 

Croucher and Clark compare all available theoretical and experimental 

( i.e . , from thermal Q data) Cg values for alkali atom interactions with 

25 

atoms and non-reactive diatomics; with a few notable exceptions the 

results are within mutual uncertainty limits. Clearly the dominant term 

in the total cross section is that due to the asymptotic Cg coefficient. 

However, a proper inversion of Q(v) data is highly desirable. Attempts 

2 6 

in this direction by several workers have not yet been put into practice. 

Irrespective of this problem (the higher-order deviations from: a 
_ ^ 

pure r -dependence in the range of r probed by thermal Q(v) measurements), 

the influence of the potential well and the short-range repulsive force is 

2 ?a 

a separate and important question. It was pointed out in 1961 that 
extrema in the velocity dependence of the cross section for atom-atom 
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scattering would thereby be expected. These "glory undulations" could be 

27c 28 29 

understood ’ from semiclassical considerations , and the extrema 

27 

indexed unambiguously. A Jef f reys-Bom (JB) (high-energy) approximation 

2 7b 

for the velocity-dependent maximum phase shift 'n was derived ' ,c 

im * 

making it possible to express the index N of the glory extremum as a 

linear function of — (valid in the high-velocity limit). Here v is 

V N N 

the velocity of the Nth glory extremum, determined from a plot of 

AQ = Q - Q vs. 1/v (where Q is the local mean of Q averaged over 

5 /2 

the undulations) or of vQ [dr (C,) ] vs. 1/v (where (C^) is 

o app ' ' 6 app 

“6 

the apparent Cg value assuming a dominance of the r term in V(r) and 

2 7 

Eq. 1 and a ^(6) from Table 1). The limiting form of (v) suggested 

-] 

N 


3 -1 

that the slope of a plot of N-g vs. v M (say I) is proportional to the 


product £r 


m 


T 



G r 

^ 1 YA 


> 


( 2 ) 


where a^ is a constant which dependd upon the functional form of V(r). 

30 

The quantity a^ was evaluated for several- realistic model potentials: 

it varied over the range 0.34 - 0.50 for the cases examined, and was 

found to depend primarily on the reduced curvature, K , of the potential 

minimum (being less sensitive to other features of potential functions). 

Thus, on the assumption of a functional form or model potential, for which 

a. was known, &x could be evaluated from the measured I or 
■L m 



7 


In an attempt to extract further information from the glory-velocities, 


Bernstein and O'Brien carried out a higher-order expansion of 


in 


30 


terms of reciprocal powers of velocity and of energy, assuming various 

model potentials. Out: of this arose an improved procedure for utilizing 

31 

precise extrema data which involved a graph (designated a BOB plot) of 

3 12 

(N-g)v vs. 1/E^ (where = ^yU-v^, whose intercept is the I of 

Eq.(3 and "limiting slope" . For the commonly used model potentials 

2 

it was found that the quantity £ r^ could be directly determined from 

S. essentially independent of tC ( to be contrasted with the relation 
i ) 

between <~r and I , via the a, of Eq.(2 ). 

The BOB procedures were criticized by several workers, however, on 

several counts. First, as recognized from the outset, itdid not take 

advantage of data on the glory amplitudes, which contain information on 
2 8 30 

the potential. Second, the deduction that a^ depended primarily 

upon the curvature K_ was empirical, based only on a few simple model 

potentials. The failure of this correlation was pointed out by D(lren 

and Schlier^ 3 by the counterexample of a flexible, multiparameter, 

potential for which K. could be varied independently of a^ . The most 

32 

penetrating analysis was that of Mason, Munn et al . who showed the 

origin of glory undulations by an optical analogy. They found that the 
extrema-spacing (and thus a^) is better-correlated with the (reduced) 
"area" of the potential well than with K , so that to a fair approxi- 
mation (i.e . , within a range of ± 20% for all model potentials tested, 
including the pathological square-well): 

oo 

I oc , (3) 
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19 b 


where cT is the (usual) zero of the potential. Third, in the case of 

the alkali-Hg systems, certain discrepancies have arisen. For Li-Hg the 

potential derived from-Q(v) data and a BOB-treatment by Rothe and Vene- 
33 34a 

klasen (as extended by Olson, whose analysis took cognizance of some 

34b 

angular distribution data ) appears to violate the empirical "Similar 

35 

Potential Hypothesis" of Stwalley, when compared with the Na-Hg and 

K-Hg potentials. A still more damaging finding was that of Buck etal. 

who compared the potential parameters of Na-Hg estimated from a BOB 

analysis of their Q(v) data with the "true values" (from a full "Buck- 

2 

inversion"), and found a discrepancy of some 207, in € r . They attri- 

m 

36 

buted this to a difference in the BOB expansion coefficients (for ^ ) 
characterizing the actual V(r.) for the alkali-Hg systems, and those for 
the model potentials of Ref. 30. At this stage it therefore becomes 
necessary to face the main question directly: given extensive data, of 
good (but finite) accuracy^ on glory extrema in the thermal energy (so- 
called "low-velocity") range, i.e. . , Q(v) as discussed, how much in- 
formation on the interatomic potential is available, i.e. , how much can 


be extracted from such data, and what is an optimal procedure to accomplish 
this goal? The present paper addresses itself to this problem. 
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II. FORMAL RELATIONSHIPS 

The total cross section Q(v) is assumed to be a sum of two components: 

a smoothly varying Q(v), dependent primarily upon the long range potential 

2 7,28 

constant C^, and AQ^(v), an oscillatory term, causing the glory extrema: 


where 


Q[v) = Q (v) -+ Afi ? M 

qm = q su .M ^ 


(4a) 


* TT 


% L, 


** H 10 


S<n ( 2. rj — 137” \ 

V I *n ^ i, 


(4b) 


where, as usual v is the relative velocity, k = yw.v/'fi is the associated 

wavenumber, L is the glory angular momentum, i.e . , the value of the 
S 

orbital quantum number Z corresponding to ^ , the maximum phase shift, 
and is the second-derivative of with respect to Z evaluated 


at Z = L 


g 


Introducing the reduced variables 


27c, d 


X' 5 l/ ir - 
A 5 4 /*■ 


r 


where r^ is, as usual, the position of the minimum of V(r) ( Eq.(4b)may 


be rewritten: 


Affl A ) • (ffr^ ^' y l3 w f - %) 

J V ' r ,2.-* i A '4 L m 'm H ) 


- *x 


A 




(5) 
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30 37 38 

Consider the following expansions ’ 3 * * based upon the high-velocity 


limit ; 


«i_ + 


+ 


5 a- -+■ 


- | , + 2l + £± (f 

\ £ ‘ £ ~ £ “ 

where E (sK) = E/£ is the usual reduced energy (£ being the depth 

of the potential well). 

30a 

It was shown that a,, B , and c. are constants determined by 

1 o 1 J 

the reduced curvature of the potential and that , 

can be taken as constants independent of the form of the 

potential and of K, £. and r , at least within a certain class of 

m 


parameterized potentials. Under the (strong) assumption that they are 

essentially "universal constants", the unknown character of , B 

2 , 2 1 (m g 

and Id /dp / . is completely specified by the values of a,. B , and C, . 

lm J /> 1 o 7 1 

The following expansion is more pertinent for the purpose of obtaining 

AQ via Eq. (5) than the expansion of p and /dfi'-Jn, of Eqs. (6b, c) 


r-^Ij r y \_ 

r* /z 

q b- i 

1 + Ml t 

tk + ... 1 

l/S. 1 * 

Jo J 

L £_* 

£**- J 


i / '/i- 

where g Q = B^/C^ . As for the A's, B's, etc.^it will be seen below 

that are essentially independent of the form of the potential, 

within the same class of parameterized potentials. Combining Eqs. (5), 

(6a) and (7) leads to a simple expression for the glory contribution: 
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(*) 

where 


cb - i-Ar •> f - — - t - (<k + + 

) - ?£ . 

(8b) 


/ H- 


The problem of evaluating Q(v) of Eq.- (4a) remains. 

write 

Formally we 

may 

S(V) + %&(*) , 


(9) 


where ^Q(v) is the "correction" for the higher order terms in the inverse 
[Jower expansion for V(r), which may be expected to be slowly-varying. 

Fully aware of the caveats mentioned in Sec. I, one can attempt to approxi- 
mate this correction assuming only the dipole-quadrupole contribution. 

13ci _g 

It has been shown previously that when the C D r term in the 

o 

potential begins to affect the cross section, its effect is to introduce 
an additive correction to Q , essentially velocity°independent , namely 

r&M = c */c<, ■ do) 

1. 6 2 39 

Values of the Cg/Cg ratio range typically from 2.7% for the He-He 

£o 2S%^ for Cs-Cs, so the Cg contribution to Q should be fractionally 

small, though not negligible (see, however. Ref. 25). As discussed in 





13 


preceding page blank not FILMED 


III. DETERMINATION OF THE EXPANSION COEFFICIENTS 

* r 2 * . 2 1 

Using the tabulated values of , (3^ , and Jd ‘ /r j m /df3 given 

37 ^ 

by O'Brien the values of the coefficients (A.B.,C.,H.) were determined 

ly I 1_ 1 

in the following way: the reference potential was taken to be of the 

standard L.-J. (12,6) form, although for the purpose of comparison the 
calculations were repeated for an Exp(^<,6) potential, with c< = 13.772. 

Using 20 tabulated points in the range 1 ^ E 100, for each of 
Eqs. (6a), (6b), (6c) and (7) a weighted (assumed constant relative 
error) forward curvilinear regression was carried out for the pertinent 
quantities ( i.e . , , etc.) with a^B^ and held fixed at 

their theoretical known values. The regression was stopped when the 
next term to be added was not significant at the 90% level (via the F- 
test). The results are given in Table II. It was found that reproduction 


to full accuracy of ^ ^ was obtained by including terms to A^ in (6a), 

L 2 * ii 

- d' Yir/dP 1 using only the term in (7). 

In terms of these truncated forms, Eq. (8) becomes 


AQ s tv)-(ZTT 0 , 


(11a) 


where 


4> s 


-3rr + Vfr, 


^7 A, ^ fSt + MX 

r n E+ J 


(Hb) 


Isolating the velocity dependence and assuming H^,A^,...A^ to be constants, 
Eq. (4) becomes (in the limit ^Q(v) -> 0) : 


Q(v)= ^ + Y x v ‘ /z + ^ P J ') *** ^ 


(12a) 
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with 



(12b) 


where 


*> * 

Y l 



fjQ (*/*)* 

A 


(13a) 

(13b) 

(13c) 

(13d) 


and € are free parameters. 

Thus, given C^, g Q , a r^, and £, then •••, may be 

calculated from Eqs. (13), using p? ( 6 ) = 8.083 from Table 1, and the 

b LL 

total cross section as a function of velocity given by Eqs. (12). Table 

III shows a comparison of Q(v) calculated from Eq. (12) for a specified 

JWKB 

L.-J. (12,6) potential with computed JWKB cross sections Q . (v)'. The agree- 
ment is exact within the accuracy of the Q^^^(v). 

Several remarks regarding the use of EqS, (12) are in order. Firstly, 
it is necessary that Eqs. (4), together with (9) for $= 0 , will 

represent the actual cross sections over the energy range of interest, 

“6 

i . e . ,only the r long-range term of the potential is significant. 

Secondly, the use of Eq. (8) requires that the expansions of Eqs. (6) 
and (7) converge sufficiently rapidly over the energy range of interest. 



(Note that the expansions are singular at E = 0 .) The coefficients 

were found to be independent of the form of the potential 

only for a specific class of two-parameter potentials and there is no 

theoretical reason to believe that they are "unique" or "universal" to 

all realistic interatomic potentials. In the general case, the best 

values of Q(v) will be obtained using and r^ instead of 'the 

true ' values £,r , where and r^ are the potential constants for 

’ nr m 

the L.-J. (12,6) potential which best fits the actual potential over 
the region of r which most strongly affects the Q(v) . 



Table II Expansion Coefficients 
Coefficient of terms in the ~nth power of E fo 
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Theivalues indicate 957, confidence interval half widths determined from fit. 

Potential A is the L.-J. (12,6) while B is the Exp( c<,6) with C<-= 13.77£y both corresponding to the same 4t 
Coefficient of variation of the fit expressed in 7,. 

Fitting was performed with this coefficient constrained to this theoretical value, from Ref. 37. 
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Table III. Comparison of approximate Q(v) from Ec 


calculation 



1.585 0.9159 1124 

1.63 0.9688 1145 

1.73 1.079 1114 

1.83 1.221 1017 

1.94 1.372 1021 

2.05 1.532 1050 

2.113 1.629 1031 

3.548 4.590 805 

5.012 9.159 681 

10.59 40.91 551 



1121 

1145 

1114 

1016 

1020 

1049 

1030 

805 

681 

552 


L.-J. (12,6) potential for which, in general, a^ = 0.421559, g Q = 0;.T86299. 

Here the specific parameters of the model are £= 13.3 cpe , r = 4.92 $ 

m ’ 

C 6 = ^ ^ r m > ~ 75* x 10^ cpe - and yU = 9,69925 x 10 ^ g [note: 

-14 

1 cpe = 1 centipicoerg = 10 erg.] 

Calculated from Eq. (12). 
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IV. FITTING OF EXPERIMENTAL DATA . 

From Eqs . (4)) (8) and (9) it is possible to construct a hierarchy of 
models for experimental total cross section data, which give successively 
more information about certain potential parameters, but at the expense of 
an increasing number of approximations concerning the form of the potential. 

From Eq. (8^ Q(v) goes through maxima and minima when sin 4 = + 1. 
These extrema are usually indexed^' 3,c in the order they occur, i.e., 

N = 1 for the first maximum, N = 1.5 the first minimum, etc., It has been 
shown^b-’ C that this requires (cf. Eq. (8b)) 

(W- i)v N 

The first two models to be presented are based upon Eq. (14) and the 

3 

"experimental" variable Y s (N - g) v N • The third model is based upon 
Eq. (12). 

A. Model I . 

It is assumed the extrema indices N and extrema velocities v have 

N 

been determined experimentally. The model for the experiment is 

Y - i + s, x x 2- -+ , d5) 

where X = E A curvilinear forward regression of the experimental data 
is to be performed, keeping all terms that are significant (at, e.g . , the 
90% level). In general, the number of terms required to fit Eq. (15) will 
depend upon the experimental energy range. 


= 




( 


a, + 


/\jJ 

E 


-f 


)• 




(14) 
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Assuming that the expansion 
experimental range so that I, S^, 
the following ansatz can be made: 



30b 


) of Eq. (6a) holds over the 
are "experimental constants" 


X - 

S, = fee 1 -r m /T X fc) A, 

> , 

s. = 6^ r„/TT*) Aj 


(16a) 

(16b) 

(16c) 


With no other assumptions, 


€ r 


m 


can be determined from: 


r m 


= US I 


( 2 ) 


Under the additiona 1 assumption that and A^ for the unknown potential 
are the same as for the "standard" potential (say, the L.-J.(12,6) form), 
then 

r-L r = ni 5_ 

C 5 (17a) 

6 = ^| S 2./^Z S / > (17b) 

which are the BOB equations (Ref. 30b) . 

B. Model II 

Suppose N and are experimentally determined as above. If, 

in addition to the assumptions of Model I, it is further assumed that 
A^,...A^ are the same as those of the particular L.-J.(12,6) potential 
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w 


ith the same £ and r , the model is 


m 


Y - I -f s; Z, + e s; Z 2 - £ 2 5,'z j + e y S,'z. 


(18) 


where Z; = A; /E J • Then the following ansatz can be made: 
J J * 


a,G r 


Ml 


ir£ f 


e z r = 

1 n\ 



(2) 


(19) 


If it is assumed that A^, . . .A^ are known ("universal constants") it is 
more efficient to use Model II than Model 1, since the lower number of 
free parameters and correct functionality of Eq. (18) leads to smaller 
confidence intervals of the coefficients 1, S^, and . It should 

be noted that Models I and II are the same for energies sufficiently 
high that the last two terms in Eq. (18) can be neglected. 

C. Model III 

Suppose that Eqs,(8) hold, sufficiently well so that and 
Ap . . .A^ are constants over the experimental energy range, namely those 
for a L.-J.(12,6) potential, and the only contribution to the non- 
oscillatory cross section is Q . Then the following model for Q(v) 

o LL 

results : 


C (v ) " Y> H- V Z ( l ■+ £ ) £-tVl $ 

4 r + Y H (w t + ety^ + 

r- v s > 


(20a) 


where 


(20b) 
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C 


where W = H, /E, and W.- = A.- /E^ for \ = 1,2, 3,4 . If the 

o i j j i 

assumptions 

are valid, then 



(21a) 


(21b) 

r 3 = 

(21c) 

Yf -- J r„/4: . 

(2 Id) 

By determining )f^, and £ from a nonlinear regression of 

Eq. (21), estimates of C, , a,, r , and g may be obtained: 

n ’ 6 7 1 nr o 



(22a) 

«, = <r 3 e/a; 

(22b) 


(22c) 

^ V J > V y / \ 3/2 - 

(22d) 


Note that if there is a multiplicative bias in the apparent 

O(v), i.e„, Q (u) = (1 + -f ) O (v) for some constant -f , then the 
’ 7 true app 

values of and Y, will be biased by exactly the same factor (1 if -f ) . 

If a small component due to the C g term contributes to Q(v), this 
may be included in the model by modifying Eq. (20a) to 


Q.(v) - ^ V ^ ^ v * ( I ^ ^ 


'i 


(23) 


where is the F(Q) of Eq. (10). 

D. Remarks 

It may be expected that the value of Model I will decrease as the 
data extends into the low-velocity range, since the confidence intervals 
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of the coefficients will become large. This reflects the deterioration 
of the expansion of Eq. (6a). 

While Model II extends the viable range of Eq. (6a) further into 
the low-velocity region, this is done at the expense of assuming that 
the potential "resembles" the L.-J.(12,6) functional form. Furthermore, 
the nonlinearity in £ of Eq. (18) will introduce a bias in the con- 
fidence intervals calculated by this model. 

Model III is the most powerful of the models in that Ch and g 

6 °o 

may also be determined (and perhaps C 0 ); it uses the full set of 

O 

observed Q(v), and requires no prior graphical analysis. On the other 
hand, the estimates are more sensitive to contributions from terms of 
the type -C /r for n>6, and the model is strongly nonlinear in several 
parameters ( £) , possibly leading to convergence problems in 

the fitting procedure, and requiring fairly good initial estimates of 
the parameters. 

In any of the. models, the use, of the tabulated values of ' 
and A^, . . .A^ for a L.-J.(12,6) potential is tantamount to fitting the 
parameters dependent upon these quantities to a "standard" functional 
form of potential. For example, the £ obtained from a fit of Model II 
will not be the best estimate of the actual well-depth of the potential, 
but the which leads to the L.-J.(12,6) potential which best approxi- 
mates the true potential. For this reason there may be differences for 
potentials which cannot be well-represented by the "standard" form. 
Furthermore, unless the true potential is close to the standard form, 
the values of £■ and r^ will be "experiment-dependent", in the sense 
that the "best" L.-J.(12,6) potential describing the total cross sections 
will be different from the best for, say, the differential cross section, 
ds- found frequently in the literature. This is, howevei; a defect of all 
calculations involving a model potential form. 
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V. Example Calculations 

Three test cases were used to provide examples of the usefulness 

of Models I and II. Case A is the set of 13 extrema data for the Na-Hg 

19b 

system of Buck et al . Cases B and C are hypothetical, representing 

o 

a L.-J. (12,6) potential with 6=8 cpe, r^ = 4 A, where + 1/2% and + 1% 

normally distributed errors have been added, respectively to the Q's. 

Test calculations were carried out for this potential with no added error, 

giving the expected results. A summary of the calculations appears in 

Table IV. Cases B and C were intended to mimic case A; for each of the 

ft 

cases there were 13 data in the range 1 < E <10. (All of the data are 

in the fairly low velocity range, where the results are expected to be 

less precise . ) It should be noted that the several estimates of a n r 

i m 

obtained for case A are of significantly better precision than the value 

obtained from the "Buck- inverted" potential for this system. 

Examples of the application of Method III are given in Table V. 

Both cases D and E were L.-J. (12,6) potentials with £r = 13.3 cpe and 
o >V 

r = 4.92 A, with 48 data in the range 1 < E < 40, with + 1/2% (relative) 
m — — — 

g 

error in Q(v) added. For both cases D and E, C, = 2 <£- r , and for case 

b m 

°2 °2 

E, 3 = 1 A ( i . e . , s ^ ~ 2 A ). From the results it is seen that, at 
least for this example. Method III is indeed sensitive to the presence 
of Cg (even when 8 is small) . 

Examination of the results (Tables IV and V) leads to the following 

conclusions: (1) the error intervals in the fitted parameters decrease 

successively from Method I through III; (2) a small amount of experimental 

2 

error induces successively larger errors in £ r^ and in <3 , when the 

ft 

data are mainly in the low E range; (3) the true values of the parameters 
usually fall within the fitted error intervals. 

Listings of FORTRAN IV programs which perform the fitting calculations 
for Models I, II and III are given in Ref. 40. 



Table IV. Example calculations of Models I and II 
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3 b 

Table V. Example of the Use of Method III ’ 





CASE D 

CASE 

E 

QUANTITY 

TRUE 

VALUE 

INITIAL 0 

VALUE 

V ° d 

v ° e 

v ° d 

C 8 * 0 e 

C.V. (%) f 

0.25 

- - 

0.26 

0.26 

0.39 

0.26 

0^x10 ^(cpe-$6) 

3.773 

3.78 

3.777 
+ .007 

3.753 
+ .036 

3.776 
+ .011 

3.752 
+ .036 

8 o 

.186299 

. 190 

.212 
+ .063 

.211 
+ .072 

.213 
+ .079 

.212 

+.072 

G. (cpe) 

13.3 

14.0 

14.0 

+1.0 

14.0 

+1.0 

14.0 

+1.5 

14.0 

+1.0 

a l 

.421559 

0.44 

.435 
+ .036 

.435 
+ .049 

.435 
+ .022 

.435 
+ .049 


4.92 

4.7 

4.5 
+ .7 

4.5 

+ o 8 

4.5 

+.7 

4.5 
+ . 8 

a l €r m (cpe "^ 

27.585 

28.9 

27.62 
+ .08 

2 7.62 
+ .08 

27.62 
+ . 13 

27.63 
+ .085 

£ 2 r (cpe 2 -&) 
m 

870 

921 

888 
+35 • 

889 

+35 

889 

+52 

889 

+35 

CgXlO 5 (cpe-^) 

0 8 

0 

- - 

3 . 8+5 . 6 

- - 

7.4+5. 7 


Fit performed with assumed constant relative error. 

+ indicate 957. confidence interval halfwidths from fit. 

Fit for Cg ^ 0 was performed using Cg = 0 parameters as initial values. 

-8 

Fit assuming no r term is present. 

“8 

Fit assuming a r term is present. 

The coefficient of variation of the fit in %. 

°2 

For case Ej C - ' *“ . 


g 



26 


VI. Concluding Remarks 

41 

Total cross sections Q(v) for the scattering of (S-state) atoms 

42 

in the low-velocity region contain information primarily on the long- 
range tail of the interatomic ( 2 -state ) potential. However, it appears 
that over the region of r probed by the Q measurements the higher- 
order terms in the inverse-power expansion of V(ry) ( i . e . , those beyond 
the leading r ^ term) contribute significantly to the potential. Thus 
it is not in general possible to extract from the smooth Q(v) measure- 
ments ( via the SLL approximation of Eq. (1)) accurate values without 

-8 

taking cognizance of the C D r and higher terms. The effect of this 

O 

deviation from asymptotic r ^ dependence over the range of r affecting 

the experiments is to yield a residual velocity dependence in the 

"apparent" calculated via Eq. (1) from the smoothed Q(v). Assuming 

that, for values of r ^R(v), where R(v) = (Q(v)/2'TT) , V(r) can 

be well-approximated by V(r) ~ — Q r 6 (| -+- /Sr" 2 -) , and provided that 

2 -8 

over this range p/r « 1 , the influence of the r term is simply 

/ O \ 

to introduce a constant additive correction § Q to the SLL-approxi- 


mated smoothed cross section. 


The glory extrema-ve locities and amplitudes are governed largely by 
the characteristics of the potential well. The results show greatest 
sensitivity to the "area of the well", which in turn is roughly pro- 
portional to the product £r^ and is strongly affected by the curvature 
K. of the potential. The extrema-spacings ( i . e . , dN/d(l/v^), in a well- 

defined limit) yield the quantity I«ca^£ r^ with good accuracy. The 

2 

primitive BOB analysis (intended to yield 6 r^ directly from the glory- 
velocity data) suffers from the disadvantage that it is tied closely to 
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the "similar potential" hypothesis, i . e . , that the wells are all some- 
what "similar" in shape, and, moreover, are similar to those of several 
commonly used "realistic" interatomic potential functions, for which 
the expansion coefficients (A^, A 2 , etc.) are essentially invariant 
("universal") constants. 

In the present paper improved procedures are developed for the 

analysis of Q(v) data. The use of the given values of the expansion 

coefficients A^, . . .A^, H^, a^, is, however, equivalent to fitting 

the data to the best "model" potential (of one of the standard forms, 

such as L.-J.(12,6), Exp(od,6), etc.). Applying the methodology to 

19b 

the data of Buck et a 1 . on N(v^) for the well-characterized Na-Hg 
system indicates that the resulting potential can reproduce the v 

N 

data to within experimental accuracy (here 0.3%). Since the thus- 
derived potential differs significantly from the best V(r) determined 
by the Buck inversion procedure employing all scattering data (including 
differential cross sections, rainbows, etc.), the total cross section 
data simply do not contain the desired information. At the very least, 
if both AQ(v) and Q(v) data can be reproduced within experimental 
error by the present procedures, then the data have been essentially fully 
exploited ; the derived potential is by definition "indistinguishable" 
from the true potential until further data (either extended energy range 
or angular distributions) can be used to remove the ambiguity. 

The literature on neutral atom-atom and atom-molecule scattering 
from the viewpoint of intermo lecul ar force determinations is extensive 
and has been reviewed by Bernstein and Muckerman in 1967 4 and by Schlier 
in 1969 . In the Appendix is presented an annotated bibliography of 

recent experimental papers which report total cross sections for systems 



involving only atoms and/or hydrogen molecules. The subject of dif- 
ferential cross sections is beyond the scope of this paper but is 

clearly an extremely important aspect of the general problem of the 

45 

experimental determination of interatomic potentials. 


ACKNOWLEDGMENT 


The authors gratefully acknowledge valuable correspondence and 


discussions with Dr. U. Buck on the subject of glory-extrema analysis. 



29 


APPENDIX 

ANNOTATED BIBLIOGRAPHY ON TOTAL CROSS SECTIONS (Since Ref. 43) 

(Systems involving only atoms with atoms or hydrogen mo'lecules) 

J. Politiek, J. Los, J. J. Schipper and A. P. Baede, Entropie 18 , 82 

(1967). K-He, Ne. High-velocity, transition region, Q(v) . 

G. E. Moore, S. Datz and F.Van der Valk, J. Chem. Phys . 46, 2012 (1967). 

q / 

He - H 2 ,D 2 ; He - He, H 2 . High-velocity, Q(v) . 

H. L. Kramer and P. R. LeBreton, J. Chem. Phys. 47, 3367 (1967). Na - Ar, 

Kr . Glories , Q(v) . 
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